Entanglement swapping between multi-qudit systems 
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We generalize the entanglement swapping scheme originally proposed for two pairs of qubits to an 
arbitrary number q of systems composed from an arbitrary number rrij of qudits. Each of the system 
is supposed to be prepared in a maximally entangled state of mj qudits, while different systems 
are not correlated at all. We show that when a set X]*=i a ò partides (from each of the q systems 
aj partides are measured) is subjected to a generalized Bell-type measurement, the resulting set of 
X]j=i( m í ~ a j) partides will collapse into a maximally entangled state. 



I. INTRODUCTION 

Recently quantum entanglement has becn recognized 
as an important resource for quantum information pro- 
cessing. In particular, quantum computation f 1 1. 12| , quan- 
tum teleportation ||, quantum dense coding |4j], certain 
types of quantum key distribution ^ and quantum se- 
cret sharing protocols Q are rooted in the existence of 
quantum entanglement. 

In spite of all the progress in the understanding of 
the nature of quantum entanglement there are still open 
qüestions which have to be answered. In particular, it 
is not clear yet how to uniquely quantify the degree of 
entanglement [0, ||, [| [ïo|, Q, or how to specify the in- 
separability conditions for bi-partite multi-level systems 
(qudits) |Q. A further problem which waits for a thor- 
ough illumination is the multiparticle entanglement 
There are several aspeets of quantum multiparticle corre- 
lations, for instance the investigation of intrinsic n-party 
entanglement (i.e. generalizations of the GHZ state [flijl ). 
Another aspect of the multiparticle entanglement is that 
in contrast to classical correlation it cannot freely be 
shared among many objeets |ïl| 0, [Ï8|, |Ï9| . 

In this paper we want to concentrate our attention on 
entanglement swapping. This is a method designed to 
entangle partides which have never interacted. The en- 
tanglement swapping has been proposed by Zukowski et 
al. |2(J for two pairs of entangled qubits in one of the Bell 
states. Zeilinger et al. [^ï] have generalized the entan- 
glement swapping to multiparticle systems. Bosc et al. 
p2[ proposed a different version of multiparticle entan- 
glement swapping and suggested a few interesting ways 
of using this phenomenon. Bose et al. investigated 
the purification protocol via entanglement swapping with 
non-maximally entangled states. This approach has been 
further improved by Shi et al. (24), and Hardy et al. 
p5[ . Dclayed choice entanglement swapping has been 
proposed and analyzed by Peres |26] . In J23, Bfj the idea 
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of entanglement swapping has been generalized to con- 
tinuous variables. The use of entanglement swapping for 
purification in continuous dimension has been proposed 
by Parker et al. [^7|. Entanglement swapping has been 
used not only for purification but also for cryptographic 
protocols (see, for instance, 28 ). Finally, we note that 



entanglement swapping has been performed experimen- 
tally by Zeilinger et al. 

In this paper we will unify all theoretical approaches 
to the entanglement swapping in one generalized scheme. 
We present entanglement swapping for systems consist- 
ing of any number of entangled systems, each composed 
of an arbitrary number of qudits (i.e. quantum partides 
with Hilbert spaces of an arbitrary dimension D). This 
new unified approach allows us to discuss in detail various 
scenarios of multiparticle entanglement. Moreover, our 
formalism applies to all possible situations when quan- 
tum systems are maximally entangled. We do not discuss 
in this paper entanglement swapping between partially 
entangled systems. 

In section II we present a relevant formalism for 
a description of kinematics of quantum states in D- 
dimensional Hilbert spaces. Section HI serves as a simple 
introduction to our swapping scheme. We show how via 
a Bell-type measurement entanglement swapping can be 
realized. This idea is extended in section IV for the case 



of two entangled states, each having an arbitrary finite 
number of partides. In section^ the most general entan- 
glement scheme is presented. We summarize our results 
in section VI. 



II. ENTANGLED QUDITS 

Let the £)-dimensional Hilbert space be spanned by 
D orthogonal normalized vectors \xk), or, equivalently, 
by D vectors \pi), k,l — 0, . . . ,15—1. These bases are 
related by the discrete Fourier transform 
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Without loss of generality, we assume that these bases 
are sets of eigenvectors of two non-commuting operators, 
the 'position' x and the 'momentum' p, such that 



x\x k ) = Xk\x k ) , p\pi) = pi\pi) 



where 



Xk = L 




Pi 




(2) 



(3) 



T he leng th, L can, for example, be taken to be equal to 
\J h/ujm, where m is the mass and ui the frequency of a 
quantum 'harmònic' oscillator within a finitc dimcnsional 
Fock space (in what follows we use units such that H = 1). 

Next we introduce operators which shift (cyclically 
permute) the basis vectors |32 : 



Rx{n)\xk) — \X(k+n)moàD) 5 

Rp(m)\pi) = b(/ +m ) mo di5) , (4) 

where the sums of indices are taken modulo D. In the 
x-basis these operators can be expressed as 



(x k \R x (n)\xi) 
(x k \R p (m)\xi) 



— Ok+n,l ; 



(5) 



Moreover these operators fulfil the Weyl commutation 
relation |33, p3l 



R x {n)R p (m) = exp( 



2n \ - 
i—mnj R p (m)R x (n) , (6) 

although they do not commute; they form a representa- 
tion of an Abelian group in a ray space. We can dis- 
place a state in arbitrary order using R x (n)R p {m) or 
R p (m)R x (n), the resulting state will be the same - 
the corresponding kets will differ only by an unimpor- 
tant multiplicative factor. We see that the operators 
R x (n) and R p (m) displace states in the directions x and 
p, respectively. The product R x (n)R p (m) acts as a dis- 
placement operator in the discrete phase space (k, l) [ j36| . 
These operators can be expressed via the generators of 
translations (shifts) 



R x {n) = eiq>(-ix n p) ; 
Rp(m) = exp(ip rn x) . 



(7) 



We note that the structure of the group associated with 
the operators R x (n) and R p (m) is reminiscent of the 
group of phase-space translations (i.e., the Heisenberg 
group) in quantum mechanics p7[ . 

Let us assume a system of two qudits each described 
by a vector in a D-dimensional Hilbert space Ti. The 
tensor product of the two Hilbert spaces can be spanned 
by a set of D 2 maximally entangled two-qudit states (the 
analogue of the Bell basis for spin-^ partides) pTfl 



\i/)(m t n)) 



1 



D-l 



fe=0 



(8) 



where m,n = 0, . . . , D — 1. These states form an or- 
thonormal basis in the space TL&7Í 



(ip(k,l)\ip(m,n)) =5k, m Si ín 



with 



D-l 



\ip(m, n))(ip(m, n)\ = J(8 / 



(9) 



(10) 



In order to prové the above relations we have used the 
Standard relation ^2^=o exp[27ri(fc — k')n/D] — DSk,k' ■ 

It is interesting to note that the whole set of D 2 max- 
imally entangled states \tp(m,n)) can be generated from 
the state |^(0, 0)) by the action of local unitary opera- 
tions (shifts) of the form 



\iP(m, n)) = Rp(m) ® R x {n)\iP{0, 0)) 



(11) 



In what follows we shall simplify our notation. Be- 
cause we will work mostly in the x-basis we shall use the 
notation \xk) = \k). In addition we will use the notation 
x Qy instead of (x — y) mod D. This serves to keep the 
derivations as synoptical as possible. Using this notation 
we can write down the maximally entangled state of two 
qudits as 



1 



**|n) |ne*))i, (12) 



where parameters k and l can take vàlues between and 
D-l. 

In general, M-particle maximally entangled states can 
be written as 



|*) = \ll>(l,fa,k2,...,kM-l)) 



j D-l M-l 

^e^|n) \neh)i 



(13) 



These partides are entangled in the sense that tracing 
out any (M — 1) partides leaves the reduced density ma- 
trix of the remaining particle in a maximally mixed state 
described by the density operator j^I. 



III. TWO ENTANGLED PAIRS 

First of all we study a simple example of entanglement 
swapping between two qutrits. Suppose we have two Sys- 
tems each composed of two entangled 3-dimensional pairs 
of partides. The two systems are not correlated at all and 
the state vector describing this composite system can be 
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entangled pair 




measurement 



entangled pair 

FIG. 1: A schematical description of entanglement swapping 
between two pairs of qubits (qudits). The qubits 1 and 2 
are measured in the Bell basis. This measurement results 
in the entanglement of the qubits and 3 which have never 
interacted directly. 



expressed as 

|*) = |^(0, 0))oi« |^(0, 1)> 23 



(14) 



= _(|00)oi + 
®-^=(|02) 23 



|ll)oi + |22) 01 ) 
f |10) 23 + |21) 2 3) 



- |00> l|21>23 
|ll> l|21>23 

|22>oi|21>23 



3(100) 01 102) 2 3 + |00)oi 1 10)23 
+ |H)oi|02) 23 + |ll) 01 |10) 2 3 + 
+ |22)oi|02) 23 + |22)oi|10) 2 3 + 



Now assume we perform a projective Bell-type measure- 
ment of partides 1 and 2 in the basis ([[2]) with D = 3. If 
the measurement yields \ip[r, s))ia for some fixed r and s, 
the other two partides collapse into the state \ip(l, fc))o3- 
This result of the measurement conditionally 'selects' the 
vectors of the form 



|n) |n)i|n')2|(n'-l)mod 3) 3 



(15) 



for n — ... 2, such that n' = n — s (mod 3) and k = 
s + 1. The amplitude of the vector |n)i|(n — s) mod 3)2 
is e i27r ™"/ 3 . It must hold that 



e i27T«0/3 e i27rn'0/3 



(16) 



Since é 1 



^u' llu 



Vu, u'eZ the equation (^6|) holds for 
/ = (— r) mod 3. The previous derivations yield that the 
state of the partides and 3 collapses into the maximally 
entangled state \ip((0 — r) mod 3, (s + 1) mod 3))o3 of two 
qutrits. 



state \ip{l, fc)) i<g>|?/>(í', k'))23- When we perform the mea- 
surement according to the basis (|Ï2| ) with D = 3 we ob- 
tain the vector \vp(r, s)) 12 ■ The resulting state of partides 
and 3 is again denoted as fc))o3- I n this case we 
are looking for the vectors of the form 



\n) \nQk) 1 \n') 2 \n' Qk'} 



(17) 



such that n! = n—k—s (mod 3), which yields k = (k+s+ 
k') mod 3. The coeficient of the vector \nQk)\\nQkQs)2 
is e l27r (""' c ) r / 3 . it must hold as before (see equation (|Ï6|)) 
that 

i27rnl/3 z27r(n— k— s)l' /3 í2tt{ii — k)r/3 i2-Knl/3 í2ttx/3 

(18) 

for n = 0, 1, 2, where e l2 ^ x / 3 will be part of the phasc shift 
of the vector fc))o3- This implies the congruence 

n(l + l'-r-ï) = -kr + kl' + sl' + x (mod 3). (19) 

The case n = gives — kr + kV + sl' + x = (mod 3), 
so the x must be chosen such that this congruence is 
satisfied. For n = 1 this leads to a relation 



ï=(l + l'- r) mod 3. 



(20) 



The extension to an arbitrary finite-dimensional systems 
is straightforward. It suffices to replace all ' mod 3' by 
' mod D ' and n varies from to D — 1. In equation 
( |ï~7| ) the generalization to £>-dimensional system gives 
us n' = n — k — s (mod D). Since n varies from 
to D — 1, we have D vectors of the form (p^). There- 
fore their linear combination with appropriate coefncicnts 
gives \tp(l, k))o3 and not only a linear combination of less 

than D distinct vectors of the form e l27rZ "/ D |n)o|/i G fc)3- 
We can now summarize our results as follows. 

Theorem 1 Suppose that \^} = \tp(l, k)) i ® \ip(l', fe / ))a3 
is the tensor product of two maximally entangled pairs of 
qudits. Let assume that the partides 1 and 2 are mea- 
sured via the Bell-type measurement in the basis (|Ï2|). If 
the measurement yields the result \4>(r, s))i2, then the two 
partides and 3 collapse into the state 

\ip{{l + 1' -r) mod D, (k + k' + s) mod D)) 03 . (21) 

This is a maximally entangled state of qudits and 3, 
which have never interacted before. 



IV. ENTANGLING TWO MULTIPARTICLE 
SYSTEMS 

Measurement of two partides 



Measuring a general state 

Let us consider now a slightly more complex situation. 
We have a system of two entangled pairs in the general 



Suppose we have two uncorrelated systems of qudits. 
The first system with mi + 1 qudits is in a maximally en- 
tangled state \ip(l, fei, ... , femi)), while the second system 
with 7712 + 1 qubits is in the state k[, . . . , k' m )). 



ooo-oo 
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FIG. 2: A schematical description of entanglement swap- 
ping between two sets of entangled qudits. A single particle 
from each set is measured. This measurement results in en- 
tanglement between the rest of the partides from both of the 
systems. 



FIG. 3: The same as in figure 2 except an arbitrary number 
of partides from each set is measured. 



Measuring more than two partides 



The state vector of the composite system then reads 
h,...,k mi )) <8> \ip(l\ !«[,..., k' ma )). Now we can 
choose two arbitrary partides (one from each of the two 
systems) to be measured using the Bell-type projective 
measurement. Due to the cyclic symmetry we can as- 
sume that the 'last' particle of each of the two systems 
is measured. Suppose that in a measurement we obtain 
a state \ip(r,s)). Therefore we are looking for vectors of 
the form 

|n)|n 9 h) . . . \n 9 k rni )\n') ...\n'e KJ (22) 

such that n' = n — k mi — s + k f m (mod D). To simplify 
the following derivations we put fco = k' — 0. Let 



k' 



Now we should determine the / and therefore 

i2irnl/D i27r(n~k m , -s+k ' )V /D _ 

OC- * 

._ i27r(n-k mi )r / D i27rnl/ D i2nx / D 



(23) 



(24) 



(25) 



It follows that 

n(l +1' -r- ï) = 

= —k mi r + k mi l' + sl' — k'„ l2 l' + x (mod D). 

As before (see equations (|ï^) and (^Ü|)) for n — 1 we have 

ï= (l + l'-r) mod D. (26) 

Once we have determined l, we can choose suitablc x to 
satisfy the case n = 0. This means the following congru- 
ence is equal to zero: 

-k mi r + k mi l' + sl'-k' m / + x = (mod D). (27) 
The resulting state is \ip(l,ki,..., fc mi +m 2 -i)), where 



ki — k{ 

h — h 1 

"'i — n·i-mi 



Ak 



i < mi 
mi < i . 



In this section we have presented a technique which al- 
lows us to produce an entangled state with any number 
of partides. 



Suppose that we are measuring the last a\ partides of 
the first system and the last 02 partides of the second 
system. We again assume the Bell-type measurement 
in the basis \ tp(r, s±, . . . , s ai + a2 _i)) describing maximally 
entangled states of the a% + ai qudits. 

Analogically as in the previous examples we are looking 
for vectors of the form 



\n) ... \nQ k mi -ai+i) ■ ■ ■ \nQ k mi )® 

®\rí) . . >' e K l2 _ a2+1 ) . . .\rí k m2 ) 



(28) 



such that nl = n- k mi - ai +i -s ai +k m2 - a2+ i for a given 
result of the measurement |^) = \ip(r, Si, . . . , s ai -|_ 02 _i)). 
Let Afc = fe TOl _ ai+ i + s ai - k„ l2 - a2+ i and k' a = 0. Now 
let us determine l. It holds that 



J2K(nl-x)/D i2-K(n-k mi - ai + 1 -s ai +k' m2 _ a2 + 1 )ï/D _ 



(29) 



gi27r(n — k mi _ ai + \)r / D ^i2irnl / D 

This leads again to the relation 

ï= (l + l'-r) mod D, (30) 

so the state of the unmeasured partides is 
\ip(ï, fei, ... , fc mi+m2+ i_ ai _ a3 )), where 



fej — fej 

k — V 



Afc 



i < mi — a± + 1 
m 1 — a\ + 1 < i. 



Theorem 2 Suppose that we have two entangled systems 
with mi + 1 and m,2 + 1 partides, respectively, initially 
prepared in the state 



I*) = |V(I,ki ) ...,A mi }<8|V'(f / ,*i, 



(31) 



and suppose that we subject the last a\ partides from the 
first system and the last a-i partides from the second sys- 
tem to a joint Bell-type measurement in the basis formed 
by vectors \tp(r, Si, . . . , s ai _|_ 02 _i)). Then the vector de- 
scribing the state of the remaining m\ + m,2 + 2 — a\ — ai 
partides after the measurement is 



\ip(ï,kl, ■ 



■ ! «-mi +m,2 + l- 



-01— a 2 )) 



(32) 



ooo-oo 
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FIG. 4: The same as in figure 3 except that many initially 
uncorrelated multi-qudit systems are considered. 



where 



k'_ 



k% — k-i 
Ak 



i—mi +ai — 1 

ï=(l + ï - r) mod D 



i < mi — oi + 1 
m\ — a\ + 1 < i 



This means that the remaining partides end up in a max- 
imally entangled state. 



V. MANY MULTIPARTICLE ENTANGLED 
STATES 

In what follows we describe the most general situ- 
ation for entanglement swapping: Suppose we have q 
systems. The jth system is composed of rrij + 1 (j = 
\,...,q) partides which are in a maximally entangled 
state \tp(P , k{, . . . , k 3 m _)). The different systems are to- 
tally factorized, so the state vector of the composite sys- 
tem reads 



i|^(F,M,...,^)). 



(33) 



(We note that superscripts do not denote the power, but 
they serve as indices.) Further we assume a multiparti- 
cle Bell-type measurement. Specifically, we consider aj 
partides from jth state, Vj € 1 ... q, to be measured si- 
multaneously in the basis 



(34) 



such that 



n 3 = n 1 - k} ni _ ai+1 - s 2 a2 + kl 3 _ a3+1 (mod D) (36) 

which in general can be cxpresscd as 

n l = n 1 - k mi - ai +i - Sa7_\ + K lt -a,+i (mod D ) 

(37) 

for Vi = 2 ... q. It remains to determine l. As before we 
have 

<? 

n i2-Kn?P /D _ i2Tr(n 1 -k mi - ai+1 )r/D i2irn 1 ï/D i2irx/D 

J = l 

(38) 

which yields 




E P ( fc «U-a 1 + l + 47-\ ~ K h -a J+ l) (mod D) 



Í=2 



(39) 




The right-hand side of the congruence ( p9| ) is equal to 
(mod D) which affects only the global phase. Therefore 
we can write 



(40) 



Consequently a set of J2j( m j ~ a i + 1) unmeasured par- 
tides becomes entangled due to the Bell-type measure- 
ment performed on the ■ aj partides. The state of the 
unmeasured partides is 

k\ 1 . . . , k mi _ ai+ i, k\, . . . , k^ lq _ aq )). (41) 

Together there are ÍX}j=i m j — aj + ij — l k and they 
must satisfy the condition 



k] — k] + n J — n 1 i < 



rrii — o, 



(42) 



The total number of measured partides is 52? =1 a j- After 
the measurement these partides collapse into one of the 
vectors (p4Í). Therefore we look for the vectors 



|n«e*« ) (35) 



l 1 )|n 1 efc 1 1 )...|n 1 efc 1 „ 1 _ ai+1 ) 

®···®|n')...|n«efeVa,+x) 



Theorem 3 Suppose we have q entangled systems each 
composed ofrrij + 1 partides (j = 1, . . . , q). Let the whole 
system is initially in the state ( |33"| ) . Let us subject the last 
üj partides from j-th (ij = 1 . . . q) system to the Bell- 
type measurement in the basis formed by vectors (|34|). 
Given the result of the measurement (^4|) the ^2j( m j — 



6 



dj + 1) unmeasured partides collapse into the maximally 
entangled state 



\tp(l, &x, • • • ) ^mi-ai + l; ^1; ■ ■ ■ ) Kn q -a q )) i (43) 



where 




and 



k\ = k\ + n J — n ; i < mj — <ij ; 



VI. CONCLUSION 

In this paper we have presented a general formalism 
describing entanglement swapping between multi-qudit 
systems. We have shown that by performing Bell-type 
measurements one can create entangled states (with an 
arbitrary number of partides) from particlcs which have 
never interacted before. 

Even though our formalism has been developed for 
finite-dimensional Hilbert space, it can be generalized for 
continuous variables, i.e. D — > oo. In this case qudits are 
replaced by harmònic oscillators (e.g. quantized modes of 
an electromagnètic field). Formally, in the limit D — ► oo 
we can substitute a two-qudit maximally entangled state 
by a two-mode correlated state, i.e. 



'\x n )\x n - Xk) -> \ip(x,p)) (45) 



where 



\ip{x,p)) 



'2-k 



\x)o\x ~ x)i. 



(46) 



Analogously, a multi-mode entangled state in the contin- 
uous limit can be expressed as 



M-l 



dx e ipi \x} 



(47) 



j'=i 



Once these states are defined one can formally perform 
the same manipulations as in the case of qudits, i.e. gen- 
eralized Bell measurements, etc. Nevertheless, we re- 
mind ourselves that the maximally correlated states J45| ) 
as well as (fÏ7|) require infinite energy for their creation. 



For this reason it is desirable to consider two-mode (and 
multi-mode) squeezed states which in the limit of infinite 
squeezing are equal to ( (45| ) and Jl7|), respectively. It is 
convenient to describe these two mode state in term of 
their Wigner functions. In particular, the Wigner func- 
tion corresponding to a regularized version of the state 
1^(0,0)) is f§ 

W(x 1 ,p 1 - 1 x 2 ,p 2 ) = exp i — — [(xi - x 2 ) 2 + (pi + P2) 2 } 



x exp 



,-2í 



[{xi+x 2 f + {p l -p 2 ) 2 } 
(48) 



This is a Wigner function describing a two-mode 
squeezed vacuum. If we trace over one of the modes, i.e., 
if we perform an integration over the parameters x 2 and 
p 2 we obtain from (|4q) a Wigner function of a thermal 
field where n = sinh £ is the mean excitation number 
in the two-mode squeezed vacuum under consideration. 
We note that the thermal state is a maximally mixed 
state (i.e. the state with the highest value of the von 
Neumann entropy) for a given mean excitation number. 
This means that the pure state (Q) is the most entangled 
state for a given mean excitation number. From this it 
follows that to create a truly maximally entangled state, 
i.e. the state (48) in the limit £ — > oo, an infinite number 
of quanta is needed and so infinite energy. 

From ( fÏ8| ) one can easily find the Wigner functions of 
other states \ip(x,p)). We remind ourselves that Wigner 
functions are invariant under canonical transformations 
(0). Taking into account that states \ip(x,p)) can be ob- 
tained from |^(0, 0)) by a canonical transformation (see 

&) 



\rP(x,p))=R p (p)®R x (x)\iP(0 ) 0)) 



(49) 



its Wigner function can be obtained via a simple substi- 
tution of variables from the Wigner function (^). The 
generalized Bell measurement in this representation cor- 
responds to a POVM measurement of the Artur-Kelly 
type [^6| . This formalism in the infinite squeezing then 
leads to a perfect entanglement swapping between har- 
mònic oscillators. 
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